Two component interacting bosons: 1-d exact results 
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Motivated by the experiment of two-component Bose- Einstein condensates produced in magnet- 
ically trapped S7 Rb, we study one dimensional Boson systems with repulsive <5-function interaction 
in the presence of SU(2) intrinsic degree of freedom by means of coordinate Bethe ansatz. The 
ground state and low-lying excitations are solved by both numerical and analytical methods. It is 
shown that the ground state is an isospin-ferromagnetic state and the excitations are composed of 
three elementary particles: holons, antiholons and isospinons. The isospinon is a triplet coupled to 
the "ferromagnetic" background antiparallelly. 
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Exactly savable models |l|-|l6[ are known to play an 
important role in the investigation of one-dimensional in- 
teracting many-particle systems. These approaches have 
served as a great source of inspiration for the new under- 
standings of non-perturbative phenomena of correlated 
electronic systems, for example, the spinon was explic- 
itly demonstrated on the basis of one-dimensional exact 
solution of Hubbard model 0. Among those models, 
an earlier typical one is the one-dimensional bosons with 
repulsive (5-function interaction which was solved H by 
means of Bethe-ansatz method. The method was devel- 
oped to solve the problem of spin- 1/2 fermions ||,|9| with 
5-function interaction. Since then had been there various 
extensions, such as the consideration of electrons on crys- 
talline lattice H , the generalization to the systems with 
higher symmetries 10 -12 and applications to different 
boundary conditions fl3| |lq|. Nevertheless, the features 
of SU(2) bosons with (5-function interaction, as far as 
we are aware, has not been studied hitherto. Because 
particles with "spin-1/2" are regarded as fermions, the 
two-component model with anti-symmetric permutation 
was studied thoroughly but the symmetric case did not 
absorb much attention. Recently, however, a two com- 
ponent Bose gas was produced in magnetically trapped 
87 Rb by rotating the two hyperfine states into each other 
with the help of slightly detuned Robi oscillation field 
[0,[L8[. It was noticed |19] that the ground state for 
a bosonic system can be surprisingly different from the 
traditional scalar Bose system once acquiring intrinsic 
degree of freedom where variational method was used. 
It is therefore very important to study the features of 
the ground state and low-lying excitations of the multi- 
component interacting Bose systems. 

In order to attain non-perturbative insight on the fea- 
tures of such system, in present letter we study the 
model of one-dimensional SU(2) bosons with repulsive 6- 
function interaction which is an integrable model. Show- 
ing the connection with the coupled Gross-Pitaevski 
equation [£0|,|l| related to Boson Einstein condensation 
|p7| , we first solve the Bethe-ansatz equation for bosons 
with SU(2) intrinsic degree of freedom. The distribu- 
tion feature of the quasi-momentum determined by the 
Bethe-ansatz equation is investigated in week coupling 
limit. The ground state property and low- lying excita- 
tions are studied by both numerical calculation and ther- 
modynamic limit approach. We fascinatedly found that 
the ground state is not an "anti-ferromagnetic" but a 
'ferromagnetic' state, the charge-isospin phase separation 
exists and the isospinon is a triplet in stead of doublet 
which is well known for spinons in spin-1/2 Fermi system. 

The two-component Bose gas produced in a recent ex- 
periment is known to satisfy the coupled Gross-Pitaevski 
equations: 



where G = -£v 2 + V(r) + c£ a |^ a | 2 , P = Ml/2. SI is 
Rabi oscillation field and V(r) the trapping potential. 

Considering the system trapped in a one dimensional 
ring of length L and introducing (\>\ = (tpi + ip2)/V2, 
4>2 = (tpi — '02)/v / 2 I we obtain from the above Gross- 
Pitaevski equation for real tt the following Hamiltonian: 

H = fdx[d x (f>* a ■ d x <t> a + ctttiMa ~ {-l) a Sl<t>* a <t>a 

where a,b — 1,2, the repeated indices denote summation 
and natural unit are adopted for simplicity. The fields 
obey bosonic commutation relations [(f)* (x), <pb{y) ] = 
2~2nez dab8(x — y — nL). The Rabi oscillation field con- 
tributes a Zeemann-type term. To avoid unnecessary 
confusions away from the conventional spins which dis- 
tinguish between fermions and bosons in the known spin- 
statistical connection theorem, we denote the generators 
of isospin SU(2) by J, correspondingly [1 + ,I~] = 2X Z . 
The intrinsic degree of freedom can be specified by the 
eigenvalues of T z or labeled by isospin up and down. 

The Bethe-ansatz equations for two-component bosons 
are obtained as follows. 
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where Ep{x) — (x+if3c)/(x—if3c). Eq.(fi]) determines the 
values of the quasi- momenta {kj} and the isospin rapidi- 
ties {Ay} for a N — 2M + 1 fold multiplet characterized 
by the total isospin X to t = (N ~ 2M)/2. Eq. (pi) is mimic 
to the Bethe ansatz equation of BBB case of Ref. p3] 
except a variant in the exponential. 

Equation (|l|) is obtained by the following procedure. 
Applying the Hamiltonian to the Hilbert space of N par- 
ticles and considering its first quantized version on the 
domain R \ {Py} where Py := {x\xi — Xj = 0} is the 
hyperplane defined by the (5-function singularity, we get 
that only the A^-dimensional Laplace operator remains 
in the Schrodinger operator. Thus A^-dimensional plane 
waves solve the Schrodinger equation. We sum up all 
the plane waves of which the wave vectors are just per- 
mutations of a definite k = (fei, fe, ■ ■ ■ , fcjv) according to 
the Bethe ansatz strategy. Integrating the Schrodinger 
equation across the hyperplanes, we obtain S(ki — kj) — 
[ki — kj — icP]/[ki — kj + ic], which connects the Bethe- 
ansatz wave functions between the regions separated by 
the hyperplanes and S := VS (V stands for the spinor 
representation of the permutation group Sn) relates the 
coefficients of different plane waves in the same region. It 
is worth to point out that the bosonic permutation sym- 
metry (instead of the antisymmetry) was imposed when 
solving the S'-matrix. Analogous to the case of spin-1/2 
fermions [pi , the periodic boundary condition leads to an 



eigen-equation for the product of the S-matrix. As the 
obtained S'-matrix satisfies the Yang-Baxter equation, 
the quantum inverse scattering method [ p3[ is applica- 
ble for diagonalization. After writing out the fundamen- 
tal commutation relations and evaluating the eigenvalues 
of the reference state |w) = | ft . . . f), one may immedi- 
ately realize the difference from that of spin- 1/2 fermions. 
For example, A(£)|u>) = UM ~ 6 - *c)/(£ - 6 + *c)|w), 
•D(6I W ) =H.i(€-€l)/(t-kl+ic)\u) in the notion of ref. 
p3Lp4|. Consequently, we got the result of eq. (fy). It 
is an facilitating picture that the Bethe ansatz strategy 
implies the existence of infinitely many constants of mo- 
tion, y~^j k'j = constant, in addition to the usual energy 



This requires that the smallest ki must be smaller than 
the smallest A„, otherwise the left hand side would be 
—A. Eq. (Eh also implies 
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^][sgn(A 7+ i - h) - sgn(A 7 - hi)] 

= 2(J 7+1 - J 7 + 1). (5) 
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Thus, for J 7 +i — J 7 = m, there must exist exactly m + 1 
solutions of ki satisfying A 7 < fc; < A 7+ i. Furthermore, 
from the first equation of eq. (0) we obtain 
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The logarithm of eq. (Qj) gives rise to 



ki = -j-lj + lJ2 ° l( - k 3 ~ k ^ + iJ2 e -i/2(% - K), 
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where Qp(x) := — 2 tan _1 (a;//3c), and both the quantum 
numbers Ij and J 7 take integer or half integer values 
depending on whether A — M is odd or even. As a com- 
parison, the Bethe-ansatz equation for spin-1/2 fermions 
is not only lack of the first summation but also has an 
opposite sign in the second summation in the first line of 
eq. (|2|). The momentum is easily obtained from eq. (|2]), 
P = J2iki = (EiIi-E v Ju)2tt/L 

It is inspirable to analysis eq. (||) in the strong and 
weak coupling regimes. For strong interaction c — > oo 
the wave function vanishes for any Xi — Xj and hence 
the bosons avoid each other like fermions which is in 
agreement with the discussion of quantum degeneracy 
in trapped ID gases ||jj. 

Due to @i(x) — > — 7rsgn(a;) and ®—\/2{%) — ► 7rsgn(x) 
for x 3> 1, eq. (0) for weak coupling limit c — > becomes, 

N M 2ir 

kj + jY s & n ( k j ~ k i) ~ J Y sgn ( fc J ~ Al/ ) = T /j ' 
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^sgn(A 7 - h) -J^sgn(A 7 - A„) = 2J 7 . (3) 
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The subscript of the isospin rapidity A 7 can be chosen in 
such a way that J 7 is ranged in an increasing order, then 
the second equation of eq. (g) turns to 



N 



Y s g n ( A 7 - h) = 2 J 7 + 27 - M - 1. 



(4) 
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Because |J 7 | < (JV - M + l)/2 for a given M and 
M < N/2 due to the restriction by Young tableau, the 
minimum value of the right hand side of eq. (Q) is — 7V+2. 
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L (I J+1 -I J -1). 

Obviously, for Ij+i — Ij = n, there will be fcj+i — kj — 
2nir/L if there is a A 7 such that kj < A 7 < /%+i, oth- 
kj = (n — 1)2tt/L. So an isospin rapidity 
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of value A M always repels the quasi-momenta away from 
that value. As a result, an existing A M will suppresses the 
density of state in fc-space at the point k = \„. The more 
isospin rapidities there are, the higher the energy will be. 
Thus the ground state of SU(2) interacting bosons is no 
more a SU(2) singlet but an isospin "ferromagnetic" state 
which differs from Fermi case greatly. 

For N particles, the ground state is characterized by a 
one-row JV-column Young tableau [N] , of which the quan- 
tum numbers are {Ij} := {-(A - l)/2, ..., (N ~ l)/2} 
and {J 7 } = empty. For this state eq. (||) reduces 
to the case of 0. Here it is a (N + 1) fold multi- 
plet with I 2 = N(N + 2)/4. The density of states 
per length of the ground state is plotted (fig|y left) for 
various couplings. The "particle" -hole (or called holon- 
antiholon) excitation is defined by the quantum num- 
bers: h = -(A - l)/2 + 5 ldl (for 1 < ji < A), I, = 
1,-1 + 1+^ (for j = 2, ..., N-l), and |Ijv| > (A+l)/2. 
Fig. (g) the left is the corresponding excitation spec- 
trum. The isospinon-holon excitation is characterized by 
the Young tableau [A — 1,1], i.e., M = 1. Compar- 
ing to that of ground state, the quantum number {Ij} 
changes from half-integer to integer or vice versa, ac- 
cordingly, h = -A/2 + S hjl (for 1 < j x < N + 1), 
Ij = Ij-! + 1 + 5 j . J1 (for j = 2, ..., JV), while J 1 =I 1 + n 
so that I\ < J\ < In . This is a A — 1 fold multiplet with 
X 2 = N(N — 2)/4. The excitation spectrum is plotted 
in Fig. (Q). Two branches of the quasi-particle excita- 
tion was recently shown for two-component condensate 
pa by means of the techniques in nonlinear optics. The 
density of state for J\ — 0, j\ — 1 is plotted in Fig. 
(|l|) . By comparison to the ground state where no isospin 
rapidity exists, a rift emerged at the position of the ex- 
isting isospin rapidity for small c that is consistent to our 
previous analysis for weak coupling limit. 



In the thermodynamics limit, the Bethe-ansatz equa- 
tions give rise to the following integral equations for 
the density of roots and holes respectively in quasi- 
momentum and isospin rapidity spaces: 

1 r Q 

p(k)+p h (k) = — + dk'pik'^k-k 1 ) 



dX'a{X')K 1/2 {k-X'), 



a(X)+a h (X) = dk'p(k')K 1/2 (X-k') 
-Q 
-B 

dAV(A')-Ki(A - A'), (7) 

-B 

where K^x) = n~ 1 pc/(p 2 c 2 + x 2 ). The Q and B 
are determined consistently by J_q p(k)dk = N/L and 

J-b CT (A)rfA = M/L. It is easy to check by Fourier 
transform that the state of B = oo and ah = is an 
isospin singlet, but further calculations show that it is 
not the ground state. The real ground state, however, 
corresponds to a — pt — in eq. (R), which concludes 
that the ground state is an isospin "ferromagnetic" state 
which agrees with the result of mean field theory p9[ . 
Two-particle case is a pedagogical example to under- 
stand it. For the two-body Schrodinger equation in the 
center-of-mass frame, the permutation of particle coor- 
dinates becomes the parity reflection of their relative co- 
ordinate. The oscillation theorem in quantum mechanics 
tells that the spatial wave function without zero nods 
which is also an even parity solution yields the lowest 
energy. If the particles possess a SU(2) intrinsic degree 
of freedom, their intrinsic wave function must be sym- 
metric (anti-symmetric) to keep the total wave function 
with the lowest energy being symmetric (ant-symmetric). 
Then we may easily understand that the ground state of 
Bose system (Fermi system) is of "ferromagnetic" (anti- 
ferromagnetic) . 

The highly degenerate ferromagnetic ground state 
which is merely true for the vanished Zeemann term will 
suddenly slit up into Zeemann sublevels once the exter- 
nal field is applied. The ground state whence becomes 
a polarized state when the Rabi field which breaks the 
SU(2) symmetry is turned on. 

In order to evaluate excitation energy we let p{k) = 
Po{k) + p\{k)/L (p refers to ground state). In the pres- 
ence of isospin degree of freedom, there will be holon- 
isospinon excitation in addition to the holon-antiholon 
excitation. The later is created by a hole inside the quasi 
Fermi sea fc G [— kp, kp] and an additional k p outside it, 



pi(fc) +S(k-k)= / dk'p x Ki{k - k') + K 1 (k- k p ). 



The excitation energy consists of two terms: AE 



J k 2 pdk + k 2 = Eh{k) + e a (k p ), the holon energy Eh and 
antiholon energy s a (k p ) = —Sh(k p ) are given by, 

/k F 
k 2 p\{k,y)dk, 
-k F 

p h 1 (k,y)+K 1 (k-y)= f " dk'K 1 (k-k')p'l(k',y). (8) 

J-k F 

Flipping one isospin corresponds to adding one isospin ra- 
pidity to the background of "ferromagnetic" ground state 
which brings about one hole in the fc-sector inevitably. 
The excitation energy AE — J k 2 p\dk is solved from 

pi(k)+5(k-k) = f F dk' K 1 {k-k') Pl {k')-K 1/2 {k-X), 

J-k F 

consequently, AE = Eh(k)+£i{X). The Sh is given by eq. 
(0) and Si by £;(A) = J k 2 p\(k, X)dk with 

p\(k,X) + K 1/2 (k-X)= I F dk'Kx{k-k')p\{k',X). 

J —k F 

In conclusion we found the existence of three ele- 
mentary quasi-particles: holon, antiholon and isospinon. 
From the asymptotic behaviors of those basic modes as 
k, k p , and A tend to kF, we find that both the holon- 
antiholon and holon-isospinon excitations are gapless. 
The related dispersions for finite N are plotted in Fig. 
(pt). Differing from the spinons in a Fermi system, the 
isospinon here is a triplet which always couples to the 
"ferromagnetic" background antiparallelly. Although it 
always accompanies charge excitation (holon) for single 
or odd number of isospinons, the isospinons can be ex- 
cited in pairs without exciting the U(l) charge mode. Be- 
cause of the coupling between charge sector and isospin 
sector in eq.(|7J), both cases brings about changes in quasi- 
momentum distribution and hence turns out excitation 
energy. The charge-isospin phase separation predicted 
by mean field theory |27| clearly occurs in present case 
due to the structure of eq. (|7|) . The holons and antiholons 
are quasi-particlcs created in the momentum space while 
the isospinon behaves like dark soliton J28| in the isospin 
sector that is tended to diminish an unit value from the 
total isospin X eigenvalue. 

Elaborating an experiment for two-component Bose 
gases whose transverse excitations are frozen out and the 
dynamics becomes essentially one-dimensional |29[ is ex- 
pected to check the afore-mentioned properties by the 
detection of the isospinon excitations and the measure- 
ment of excitation spectra. 
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FIG. 1. The density of state per length in fc-space for the 
ground state (left) and for the state in the presence of one 
isospin rapidity A by choosing J\ — (right). The distribu- 
tion changes from a histogram to a narrow peak gradually 
for the coupling from strong to week c = 10, 1, 0.1, 0.01. The 
calculation is made for L = 40, N = 40. The left panel is 
mimic to the fig. 2 of ref [7]. 



FIG. 2. The holon-antiholon excitation spectrum (left) 
and holon-isospinon excitation spectrum (right) calculated for 
L=20, N=20 and c = 10. 
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FIG. 3. The left figure is the dispersion relations of 
isospinon for different coupling constants where the curves 
from top to bottom correspond to c = 1, 10, 20, 40, 80 respec- 
tively. The right figure is of antiholon, holon and isospinon 
for c=l. They are calculated with L = 20 and N — 20. 



